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Abstract 

We present a calculation of the 7 d — > K'^K^ n p reaction with the aim of seeing if the 
experimental peak observed in the K~^n invariant mass around 1526 MeV, from where evidence 
for the existence of the Q'^ has been claimed, can be obtained without this resonance as a 
consequence of the particular dynamics of the process and the cuts applied in the experimental 
set up. We find that a combination of facts leads indeed to a peak around 1530 MeV for 
the invariant mass of K~^n without the need to invoke any new resonance around this energy. 
This, together with statistical fluctuations that we prove to be large with the statistics of the 
experiment, is likely to produce the narrower peak observed there. 

1 Introduction 

The discovery of a peak in the 7 ^^C — !■ X reaction in the K^n invariant mass spectrum, 

which was identified as a signal for a pentaquark of positive strangeness, the 6+ [1], created a 
turmoil in the Hadron Physics Community and stimulated many theoretical and experimental 
works. A large number of experiments were conducted, many of which found also the signal, 
but gradually other experiments were performed that did not produce the expected peak. The 
number of theoretical works stimulated by this finding in the following years was of the order 
of one thousand and, whether they could or could not find a justification for the state, they 
gave collectively a big push towards the understanding of the hadronic spectra. A comprehensive 
review of these developments was done after some time in P], were one can see the relevant 
literature on the subject. More recently a new experiment was done at LEPS on a deuteron 
target, and with more statistics, and a clear peak was observed around 1526 MeV in the K^n 
invariant mass distribution j3]. While quite some evidence is given for the existence of the peak 
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in [3], it is not ruled out that it could be a consequence of the peculiar dynamics of the process, 
with the experimental set up and the cuts applied to see the peak, rather than the signal of a 
genuine new resonance. The fact that the same reaction measured at CLASS with ten times 
more statistics and complete kinematics fails to see the 0"*" peak [¥J, is certainly an incentive 
to search for some alternative explanation of the peak observed at LEPS. A first thought comes 
to our mind: in the reaction of [3], 7 d — > K~^K~ n p, the main contribution comes from (f> 
production, which is removed by particular cuts done in the K~^K^ invariant mass in [3]. This, 
together with the fact that the particles are detected in the forward direction could lead to some 
unwanted structure in the mass distribution of K'^ n. It should not be overlooked the fact that 
in [5] only the K'^ and the K~ are detected, but not the p and the n. The reconstruction of the 
if"*" n invariant mass is done using a prescription to eliminate as much as possible the effects 
of Fermi motion and to have the p as a spectator. We shall show in the present paper that the 
prescription works to show peaks coming from genuine resonances, but also creates a distorted 
spectrum. On the other hand, the procedure does not work in the reverse direction, and if a 
peak appears it is not necessarily an evidence of the presence of a resonance. 

There is another issue which requires attention. The experiment is done in deuterium, since 
one is interested in having the n invariant mass. This means one necessarily will have 
rescattering of the produced kaons. Hence, assuming that a (/) is created and one of the kaons 
rescatters from a second nucleon, the rescattered kaon changes energy and direction and the 
reconstructed K~^K~ invariant mass does no longer match the original (p from where the two 
kaons came from. Since the (p production counts for the largest part of the cross section in the 
experiment of [3], it could be that this part of the amplitude accounts for a relatively sizable 
part of the total amplitude and through interference with the one body processes leads to a some 
kind of signal which might be misidentified in the analysis. 

There are, thus, some reasons that make us think that some peak could be produced, not tied 
to the existence of a resonance but to this peculiar combination of dynamics, multiple scattering 
and cuts. The present paper presents a thorough simulation of the physical process, using a 
theoretical model that accounts for (jj production and rescattering of the kaons. The production 
of the A(1520) is also considered and the procedure to filter events with a nucleon as spectator 
is tested with this resonance. The rescattering part versus the one body part of the amplitude 
is done accurately, using realistic deuteron wave functions and amplitudes for KN — > KN and 
KN — > KN constructed using chiral unitary dynamics. 

The results obtained show that a peak around 1530 MeV for the invariant mass of the K'^n 
pair appears, even with large statistics, without having any 6"*" term in the amplitude. The peak 
obtained is broad. However, with the statistics of the experiment of [3], fluctuations appear often 
which can produce a narrower peak in that region, very similar to that obtained in (3j . We also 
show that it is possible, after some strength appears in that region with a predetermined set up, 
to obtain a sharper peak by changing a bit this set up. All these findings give us strong reasons 
to think that the peak seen in [3] obeys to this combination of facts related to the experimental 
set up and is not a signal for a new resonance. 

The work proceeds as follows. In section [5] we explain the formalism to calculate the cross 
section for the process — >■ K~^K~np with detail and the model to take into account the 
and A(1520) productions. In section [3] we show the K^n invariant mass distribution obtained 
with the same set up than at LEPS . In section S] we show the results that we obtain when we 
introduce the A(1520), in section [5] the rescattering contribution and in section |6] we make a 
statistical analysis and compare with the results of LEPS, after which we draw our conclusions. 
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2 Formalism 



When trying to simulate the — )■ K^K^np reaction one must consider the basic features 
learned from the experiment in [J . The most important contribution comes from (f> production, 
the elementary reaction — >■ (j)p and 771 — >■ (/>n. We will implement this (fi production in our 
simulation through a minimal model which incorporates the basic structure of K~ production 
correlated by the propagator. We will also take into account empirical information from [S], 
from where one induces similar cross section production from proton or neutron targets. 

Another element in the dynamics of the — >■ K~np reaction is the production of the 
A(1520). Although its contribution to the process is small compared to that of </> production, we 
will take advantage of it to test the method used in [3! to eliminate Fermi motion in the reaction 
on deuteron. 

An important ingredient in our approach is the rescattering that unavoidably occurs in the 
reaction in deuteron. Indeed, consider a first step of 7p — 0p — > K~^K~p and then a rescat- 
tering of the K'^ or the K~ with the neutron. This double scattering step can be evaluated 
theoretically, and its ratio to the first step 7p — > (j)p can be calculated rather reliably. This is 
so since uncertainties in the elementary (p production model cancel in the ratio and the rest 
roughly requires the knowledge of the deuteron wave function and the scattering amplitude of 
(K^) with the nucleons, which one has rather under control from studies of the KN and 
KN interaction using unitary chiral dynamics [6Hll| . 

The consideration of the rescattering is important in principle in this simulation. Indeed, 
the cj) production takes the largest strength of the 7p K^K~p reaction in the energy region 
of LEPS. Yet, with the cuts of LEPS one removes this contribution, that would show up as 
a large background on the K'^n spectrum precluding the observation of some smaller signal. 
However, if one has 4> production and rescattering of one of the kaons, this latter kaon will 
change energy and direction upon rescattering, such that when the cut on the K~ invariant 
mass is implemented, these events are not removed. The kinematics of these events coming from 
a primary (j) production is very peculiar, such that it is not unlikely that they could lead to some 
peak in the K^n invariant mass distribution which could be misidentified as a signal of some 
resonance. 

The model that we take is depicted diagrammatically in Fig. [T] We have there production 
on the p plus the same with rescattering of the or the . Similarly we also have 4> 
production on the neutron plus the same with rescattering of the or the K~ . In the case 
of 7p — > K^K^p we also can have the K^p making the A(1520), which we will also take into 
account. The technical details follow in the next subsections. 

2.1 The impulse approximation 

In order to calculate the cross section for the process ^d — > K^K^np we need to determine the 
S matrix or T matrix for this reaction. Let us consider the first diagram in Fig. [1] corresponding 
to the impulse approximation, and let us evaluate its contribution to the S matrix. Considering 
plane waves for the wave functions of the incident photon and the particles in the final state. 
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Figure 1: (Color online) Diagrammatic representation of the model for the — t- K^K np reaction. 



the 5-niatrix for the diagram of Fig. [T^ can be expressed as 

1 



d'^xi- 



V 



(1) 



where Lpi ii — f,2) are the wave functions of the nucleons in the deuteron, p\ is the energy of 
the proton inside the deuteron, pp, p„, g, q' correspond to the four momentum of the proton, 
neutron, and of the final state, respectively, and k to the four momentum of the incident 
photon. The normalization factor Ni is given by Ni = Mi / Ei for baryons and Mi = 1/ \J2uji 
for mesons. The plane waves are normalized inside a box with a volume Y . 

The calculation of S"^"-* requires the evaluation of the integrals in the x\ and X2 variables. In 
order to perform the integration in x\ we separate the temporal and the spatial components and 
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write 

The integration over the time component gives the delta function for the energy conservation 

dx0^^ip'^,W+r~k«-p<l)x« ^ 2nSik°+p°^pl -q°- g'"). (2) 

In the impulse approximation the energy of the spectator, in this case the neutron, does not 
change, i.e., p'^ = P^, with p^ the energy of the neutron inside the deuteron (omitting the small 
binding energy). Therefore, we can write the 5 in Eq. ([2]) as 

5{k' +Pl-pl-q°- q'°) = S{k" +p'1+P°2-pI- pI - <f - q'\ (3) 
The evaluation of the spatial integrals is done by making the usual transformation 

R = + X2) 

f^Xi-X2 (4) 

with R being the center of mass position of the pn system (assuming an average mass for the 
proton and the neutron) and r the relative coordinate between xi and X2- In this way, the term 



d^xi J d^X2e-'^P''+^+^'-''^^'e-'P"^'JViJ\f2Mxi)Mx2) (5) 
present in Eq. ([1]) can be written as 

J d^Re-'^Pp+^+^'-^-P-^+P"')^ J dVe"5(Pp+9+9'-^-P")*'_^A/-^(^(r=*) (6) 
where we have introduce the deuteron wave function (p^r) for the relative motion 

AfiAf2M^i)M^2) ^ -^Ude^P^'^^ir) (7) 



with pd the momentum of the deuteron. The deuteron wave function is normalized as 

rfV|(p(f)|2 = 1. (8) 



We only consider the s-wave component of the deuteron wave function, which in the rest frame 
of the deuteron can be parametrized as [12] 

11 Q_ 

(p{f) = — -exp{~mir). (9) 
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The coefficients Ci and mi can be found in [13 . The integration over r in Eq. ([B]) results into 
the Fourier transformation of the deuteron wave function, while the integration over R generates 
the delta function of the total momentum conservation 



J dPRe-'^P^+^+^'-^'P-'+P-^^ = {2TTfS^{pp + q + q' -k-pd+ Pn). (10) 

Taking the deuteron energy as the sum of the energies of the nucleons, i.e., p° = pi + P2, we 
have, collecting the previous results. 



'1/3 

X S^iPp + Pn + q + q" ~ k - pd)v 



3/-. , - , - , r -..-I k + Pn-Pp-q-q" 



Using the momentum conservation law for the reaction 7c? — >■ K~np 

k + Pd = q + q" +Pp+Pn (12) 

the argument of the (p function in Eq. (jlip can be put as 

~(k+Pn-Pp-q-q'\ Pd\ 

\ 2 J = ^l^P" - yj ■ 

and, therefore, Eq. ([TT|) reads as 

5^"' = -l^^fd^f^^fp^fn^fK^^fK~ {2T:)H\k +Pd-Pp-Pn-q-q') 

>^^(pn-^y;U^K-p (14) 



Calling 



we arrive to 



A=^(^llM}j{27rrS\k+pd- 



Pp~Pn-q-q') (15) 



S('^^=-^AMpn~^]t[;^^^,^-p. (16) 



Similarly, it is straight forward to show that 



s^''=-^AM'f~PpY:LK.K-.- (17) 
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2.2 The rescattering terms 

Let us consider now the diagrams where a kaon rescatters with one of the nucleons of the deuteron 
as shown in Figs. [T};-[Tf. Analogously to the single scattering diagrams evaluated in the previous 
section, the S matrix contribution for the diagram Fig. [TJ; is given by 



if ' 



(18) 



Separating the temporal and the spatial part in the integration over the variables Xi and X2 , we 
have, for the temporal part. 



(19) 
(20) 



rf^Og.(,»+p»+gL-fc"-p\>\' = (27r)5(fcO - - ql,) 

Considering ^ P1+P2 and using Eqs. (fT9|). ([20l) . Eq. ([T8| can be written as 



X i 



t 



(21) 



where we have made use of Eq. (fT9| in order to determine the value for q^^ . Writing Eq. ([2T 
as a function of the R and r coordinates defined in Eq. (|4]) we have 



X t 



(c) 



17)6 y (2^^ 1 1 Pp Pn) 

1 



(22) 



with j an index running over the final p, n, K'^ , and the initial photon and deuteron. The 
integration over R generates a delta function for the momentum conservation and the integration 
over r the corresponding Fourier transformation of the deuteron wave function ^p{r). Taking into 
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account this fact, Eq. can be expressed as 



Pd 



(p\q +Pn - — - qex\t 



Ac) 



■yp~^K+ K p K K n 



(qex)'^ - {qexY ~m% + ie 



(23) 



where we have used Eq. (|12l) for the argument of (p. In terms of Eq. ([15]), Eq. ((23)) reads hke 



S'-"^ = -lA 



(27r) 



-tp\q +Pn - ~ r. 



Ac) 



(c) 



7p— K p K 71— 71 



{ql^Y-[qexY-ml + te 



(24) 



Proceeding in a similar way, it is possible to find the contribution for the remaining diagrams 
of Fig. [TJ which we just write bellow 



S^'^"> = -lA 



i2ny 



~ I -^1 ^ Pd -t 1 ,lu,: 



Ad) 



Ad) 



(25) 



5'(«) = -iA 



{2n) 



~i ^, ^ _P£ _ ^ I 

^V\q+Pn 2 'Jea; ] I^p^K + K-p'-if + n^A:+r 



l-ic) 
7P- 



Ac) 



(26) 



(2^)= 



Pd 



-<^\ -q-Pp + ^ + qex]t 



Af) 



if) 



-yn-^K+K-nK+p^K+p 



1 



(&)^ - (qcx)^ -mj^+ie 



(27) 

The total S matrix will be given by the sum of the contribution of the different diagrams 

/ 

S = l + J2Si (28) 



and the total T matrix is obtained using the relationship 

l-iAT. 



(29) 
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2.3 Basic model for production and the KN , KN amplitudes 

As shown in the former subsections, the calculation of the S matrix involves the evaluation of 
the amplitudes for the processes -> (t)N — > K^K^N, where N represents a nucleon, i.e., 
N ^ n, p, and KN — > KN, KN KN. In order to determine the amplitude for the process 
"/N (j)N — >■ K^K^ N , we are going to assume a minimal structure which satisfies gauge 
invariance. Let us first consider the process ^N — > 4)N (Fig. [5]) . 




Figure 2: Diagrammatic representation for the process 7p (pp. 

The amplitude which describes this process must be a function depending on the four mo- 
mentum of the photon, fc, and the p^, and their respective polarization vectors, and e^. 
Imposing gauge invariance, this amplitude must be of the form 

t^N^<t,N ^ (e-yC^fc • P0 - p^e^ke^). (30) 

If we work in the Coulomb gauge for the photons, p^e-y = —p^e^, and considering k-p^ ^ k^M^, 
we can neglect the second term in Eq. (pO)) as far as jp^l << M^. The particular structure of 
this amplitude is not essential in order to see the effects that the cuts of |3] can produce in the 
M^+„ invariant mass distribution. Since we will evaluate the cross section in a reduced range 
of the photon energy, fc*^, we can further assume this to be constant and thus 

t-)N^<t>N = OAre^e^ (31) 

where ajv is considered as a parameter which determine the magnitude of the amplitude for this 
process. Once the is produced, it can decay into a K^K~ pair as shown in Fig. |3l Therefore, 
naming fci and k2 the four momentum of the K'^ and K~ , respectively, coming from the (j) decay, 
we can describe the process 7iV — (fiN K^K~ N by means of the amplitude 



tjN^K+K~N = aNe'^tf^^— TTJ— -T7-f^5e0(fci - k2)u (32) 

with the width of the (j) and g the (f>K^ K^ coupling. Summing over the polarizations of the 
cj), A^, and using that 



^^^'^ = -9'' + -^ (33) 

A<4 'P 
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we have 



t^N^K+K-N = aN-r, —j-T^ i\t2 , -nj t ^1^^'^ ~ ^2)/x (34) 

where ajq = ajsig and we have used that = ki + k2 and taken an average mass for the kaons, 



i.e., ^\ — ^2 — K ' 

In the Coulomb gauge, Eq. ([M)) reads as 



t-yN^K+K- N = t-iN^K+K- N^i (35) 



with 



tjN~^K+K-N = ^(^NJ7 , , X2 S72~"^7F1^ ~ ^2) (36) 

(fcl + k2Y - + iM^i 

and this is the expression which we use in order to describe the process — > (f)N K^K^N . 
Note that the values of fci and ^2 depend on the diagram under evaluation. 

For the one body process fci, k2 correspond to external variables which are known in the 
integration over phase space. The two body contributions are relatively smaller, even when the 
(j) cut is done, as will be shown in section[51 There one of the momentum is an external momentum 
while the other one corresponds to qex, an integration variable, which we approximate in the 
matrix element by its value for = p2 = in the laboratory frame, where the deuteron wave 
function in momentum space has its largest value. In this form, the amplitude t^N^K+K- n can 
be taken out of the integrals in Eqs. ((24|) - p7|) . 

The t-matrix for the KN and KN interactions are calculated by solving the Bethe-Salpeter 
equation in a coupled channels approach, taking as potential the lowest order chiral amplitude as 
done in [S1|71[TT]. For the KN system and coupled channels, the scattering amplitude calculated 
with this approach reproduces the total cross section of K~p to all channels, and gives rise to 
two A(1405) states as dynamical generated resonances in the meson-baryon interaction |8l[9], 
which can be seen in the invariant mass distribution of ttS in different reactions [51[T3]. The t 
matrices for the processes KN — ^ KN and KN KN depend on the invariant mass of the 
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interacting particles, which are given by the foUowing expressions according to the diagrams of 
Fig. [T] considered 



s(^) = {Pn + q'f = Ml + ml + 2plq'° - 2p„g', 
= (pp + q'f = Ml + m\ + - 2ppq', 

s^'^ = [Pn + qf = Ml + ml + 2plq" - 2p„g, (37) 
sU) = + ^ Ml + m\ + 2pOgO - 2ppq. 

Since these invariant masses are defined in terms of external variables, the KN and KN 
amplitudes can be taken out of the integrals in Eqs. ((24)) - (|27p . 
In these way we can define the following function 



which depends upon q^^ and Q , such that from Eqs. (p4)) - (|27)) and Eq. pQ]) we can write the 
T matrix from the different diagrams in Fig. [T] as 

T'-' -4k'+pl- q' -P%q'+Pn - §)t^;lK.K-/K'-u-.K~J. ^ ^'^^ (39) 



T'^' - +pI - q"-pU+P, - %ii\K^K-/Icl,^K./. - ^^^^< 



-7 J 



7' 



■7- 



The total T matrix will be given by the sum of T^'^ i — a,b, . . . , f. We have then from the 
different terms in Eqs. (j39p 

T = ^e-y (40) 

where 

/ 



7^ = j]7^W. (41) 
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2.4 Evaluation of the cross section 



In order to calculate the total cross section for the process — >■ K~^K~pn we are going to work 
in the center of mass system, where we denote the variables with a tilde. In this frame, the cross 
section will be given by 



r f f r d-% 1 1 Mp Mr, 

J {2n)3j (277)3; (27r)3y (27r)3 2w(| ) 2w'(|') 4(ft>) ^-(Pn) 

X {2n)''5\k+pd-Pp-Pn-q-q')J2\T\^ (42) 



s-Mi 



with Md the deuteron mass, u}{q ) = ^q'^ + m\, Cb'{q') = {q')^ + m\, Ep{pp) = ^ (pp)^ + M' 



and Enipp) = y (p„)^ + M^. We are using an average mass for the proton and neutron, i.e., 
Mp = Mn = Mm where Mm = {Mp + Af„)/2. Since the momentum of the neutron is not 
measured at LEPS, we can use the delta of Dirac to remove the integration in this variable. In 
this way we have 



M, 



(Pq f d?q' f d^pp 



Mr. 



M„ 



MjJ {2^f J {2nY J {2nY2Co{q)2u'{^')Ep{Pp)E^{-pp 



q') 



X {2'K)5{y/s - Cj{q ) - Cj' {q ) - Ep{pp) - En{-Pp ■ 



q'))Y.\T\ 



(43) 



where we have used that in the center of mass frame k + pd = 0. The 5 function is used to 
integrate over the angle between pp and q + q' which is then given by 



cos6o = — — 



1 



{y/S - ujii ) -(v'it) - Ep{Pp)f -M%- {Ppf - (1+ i'f 



2|Pp|-|g + g' 

and the remaining integrals give now the cross section, which can be written as 



a = 



MdMnMp 



1 



!dEp [4^T.\T 

J {2^fCj{q)J {2^f[b'(q')J "^J ^\q + q'\^' 



s-Ml 4(27r)2 7 (27r)3w(|)7 (27r)3£::;'(|') 
X 6(1 - cos2^o)e(Vi - w(9 ) - - EpiPp)) 



(44) 



where the 6 function guaranties that |cos0o| ^ 1 and y/s — w{q ) — ui'{q') — Ep{pp) > 0, as it 
should be. In the Coulomb gauge, the element 



(45) 



can be written using the sum over photon polarizations 



(46) 
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as 



A 



(47) 



Although Eq. (H^ corresponds to the cross section calculated in the center of mass frame, 7^|T 



A 



is an invariant and can be calculated in any colinear frame. For this reason, for convenience, it 
is simpler to obtain it in the frame where the deuteron is at rest, i.e., pd = 0. For this purpose, 
we have to make a boost from the center of mass frame to a frame which is moving with a speed 



k+Pd 



(48) 



fco + Md fc" + Md 

with fc, variables of the photon in the laboratory frame, according to the following relation [15 

'k° + Md ^\p^■k 



Pi 



1 



Pi 



(49) 



with p an arbitrary vector in the center of mass frame and p the expression of this vector in the 
frame where the deuteron is at rest. 

The integrals in Eq. (|44|) are done using the Monte Carlo method, which is particularly 
suited to implement the experimental cuts. 

2.5 A(1520) production 

A detailed model for A(1520) production, which reproduces fairly well the data, is available 
in [16]. In this work the Drell mechanism with and K~ exchange is used to study the 
7P — >■ K~^K~p reaction, and the production of the A(1520) is taken into account by a new term 
where the resonance is excited by means of K* exchange. Another phenomenological analysis 
is done in [17] where the A(1520) resonance is assumed to be dynamically generated in coupled 
channels. Here we follow a more empirical approach, since we are interested only about the shape 
of the K~p distribution and the strength of the A(1520) production is fitted to the data. We 
assume that the A(1520) does not interfere with the ip production mechanism, and we consider 
a minimal structure for jp — ii'+A(1520) — >■ K^K~p compatible with the quantum numbers of 
the A(1520) (J^ = 3/2") and the D-wave character [17] 



with 



(50) 



(51) 



^K~p ^ -"'-'A(152U) T 2~ 

the A(1520) propagator, Ma(i52o) and rA(i520) the mass and the width, respectively, of the 
A(1520) resonance, 6a a parameter, a the Pauli matrices, and M^-p the K~p invariant mass. 
In this way 

(52) 
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Using Eq. pSj) . in the Coulomb gauge we find 



(fc ■ q)' 



(53) 



We can add a background to this amphtude of the form ^|caP, with ca being a parameter, and, 
then, substitute in Eq. pi)) 



X A V 



Since in the energy region considered the term 



(fc ■ q}' 



(fc ■ q)' 

P 



~l ^ Pd 



(54) 



(55) 



does not change too much, we can take it as a constant and absorb it in the parameter &a- This 
mean that in our case 



' + \[\bK?\DA? + \cA? 



<^ Pn 



Pd 
2 



(56) 



We shall see that the results at LEPS can be reproduced with the (j) and A(1520) production, 
and there is no need for an extra background, represented by the ca term in the former equation. 

The A(1520) is dynamically generated in the KN system and coupled channels [IT]- One 
can benefit from the study done in [17] to include the A(1520) production in the t matrix for 
the process K^p — > K^p which appears in Fig. [TJi. Following [17] 



,A(1520) _ ,7=0 



M \ 



■,m,M-m,] Y2,m-M{p) 



X C ( i, 2, 1 m'. A/ - m' j r2,™'-M(p )(-l)' 



Mtt 



(57) 



with p, p' the initial and final momentum of the , t^^^^xN ^ matrix for the KN system 
in isospin zero obtained solving the Bethe-Salpeter equations in a unitary chiral approach [HIIT], 
C Clebsch-Gordan coefficients and Y the spherical harmonics. We take p in the z direction, 
and then, the only not zero spherical harmonic is 12,0 (p)- In LEPS they measure angles in the 
forward direction. Then, taking p and p' in the forward direction 



12,0(P) -1^2,0 (P) 



(58) 



Let us consider m = m' = 1/2. Since C i 2, |; i 0, i = i/ 1, we have that 



,A(1520) _ „,/=o 



(59) 
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Using that 



\KN,I ^0) = ^\K-p + K"n) (60) 



we can determine the element t^'^^^^''^„_ as 



1=0 



,A(1520) _ ^KN^KN _ ,1=0 _ \ 
''K-p^K-p ^ 2 ^ ''KN^KN y"^) 

and close to the resonance, t^^^^xN '^^-'^ written like 

where g^iv is the coupling of the A(1520) to the KN channel, which according to is -0.54. 
This term must be added to the s-wave K^p — > K^p amplitude which appears in Eqs. ([25|). We 
mention already here that including this term to the rescattering part induces corrections which 
provide a small contribution, thus, has negligible effects in the final results. 



2.6 Set up at LEPS 

The LEPS detector is a forward magnetic spectrometer, therefore, it can only detect particles in 
the forward direction. The angular coverage of the spectrometer is approximately ± 20 and ±10 
degrees in the horizontal and the vertical directions, respectively. In our simulation we impose 
that the angle of the kaons in the final state with respect the incident photon is not bigger than 
20 degrees. 

The nucleons are not detected at LEPS, therefore, some prescription is required in order 
to estimate the momentum of the p and n in the reaction — > K^K^np and determine 
the invariant mass of K~p or K'^n. This is done using the minimum momentum spectator 
approximation (MMSA). The aim is to guarantee that a chosen nucleon is the closest possible 
to be a spectator in the reaction. The chosen spectator nucleon is the proton if one evaluates 
the K^n invariant mass, or the neutron when the K~'p invariant mass is evaluated. For this 
purpose one defines the magnitude 

Ppn = Pmiss =P-i +Pd- Pk+ - PK- (63) 

which corresponds to the four momentum of the outgoing pn pair. In the pn center of mass 
frame, the momentum of a nucleon is given by 

\y\Ml^^Ml,Ml) 

Boosting back this momentum to the laboratory frame, we will have a minimum modulus for 
the momentum of the spectator nucleon when the momentum pcM for this nucleon goes in the 
direction opposite to Pmiss- Thus, the minimum momentum, Pmim is given by 

Pmin = - \PCm\ ■ -TT ^ EcM ' — TT (65) 

IVlpn -^^J^pn 
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where Ecm — \/\pcm\'^ + M'^ is the energy of the nucleon in the CM frame. In this case, 
the momentum of the other nucleon wih be in the direction of the missing momentum with a 
magnitude 

Pres — IPmiss \ Pmin 

(66) 

In [3] the M^+n invariant mass for the reaction K^K~np is evaluated assuming the 

proton to have a momentum Pmin- Consequently, in this prescription, the momentum of the 
neutron in the final state will be 

^ Pmiss /r"7^ 

Pn = Pres ' p T (67) 

IPmiss I 

which is used to calculate the Mx+n invariant mass for the reaction 7c? — > K^K^iip in [3]. 
At LEPS only events which satisfy |pmi„| < 100 MeV are considered. This condition is also 
implemented in our simulation of the process. 

In order to remove the contribution from the cj) production at LEPS one considers events 
which satisfy that the invariant mass of the K^K^ pair is bigger than 1030 MeV and bigger 
than the value obtained from the following expression 

1020 MeV + 0.09 x [E^^ f (MeV) - 2000 MeV) (68) 

where E^f-f" is defined as the effective photon energy 

E'Jf = "^+^-»'^" (69) 



with Sx+K-n the square of the total center of mass energy for the K'^K~n system calculated 
using the MMSA approximation to determine the momentum of the neutron assuming the proton 
as spectator. In [3] only events for which 2000 MeV < E'^f^ < 2500 MeV are considered, a 
condition which is also incorporated in our simulation. This is the cut implemented when the 
invariant K~^n mass is reconstructed. When the K~p mass is reconstructed, the cut is the same 
but E^-f^ is now evaluated with the MMSA prescription assuming the neutron as spectator. 
These energies would correspond to the photon in a jn (^p) reaction with the original p (n) at 
rest. 



3 Results with production 

We consider first the production of ip on proton and neutron followed by rescattering and we 
look at the distribution of da/dM^+n which is obtained taking for the neutron the momentum 
provided by the MMSA prescription described in the former section. 

We have made the integrals in Eq. (H^ using the Monte Carlo integration method with 
different seeds and with a statistics of 20 million points in the variables of Eq. (pij) . such that 
we get about 66500 "good" points in the integral over the phase space which pass all the tests 
of the cuts. In order to take into account that the photon energies in the LEPS experiment run 
from 2 GeV to 2.4 GeV, we have generated 20 random energies in that range. 

In Fig. |3]we show the first distribution obtained with the cuts of [3] and using the MMSA 
prescription for the nucleon momenta as done in [J;. What we see in the figure is a pronounced 
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peak around an invariant K^n mass of 1530 MeV, very similar to the peak observed in [3], which 
was identified as a possible signature of the 0+ pentaquark. However, the peak has appeared 
here as a consequence of the cuts and not as a signal of a resonance which is not present in our 
theoretical model. 




1450 1500 1550 1600 1650 1700 1750 1800 1850 1900 



MK+n(MeV) 

Figure 4: (Color online) M^+„ invariant mass distribution with about 66500 "good" points. 

One may rightly ask oneself whether the peak obtained in Fig. |4]is a necessary peak coming 
from the LEPS cuts or it is a statistical fluctuation. The answer comes from running with much 
more statistics, which we show below. However, before showing these results let us mention that 
statistically we get a clear peak around this energy in about 1/3 of the runs with different seeds, 
but in all of them there is some trace of extra strength in this region. 

Since one may argue that statistically it is more probable not to find a clear peak (about 2/3 
probability), we find opportune to make the following test. Let us imagine that we had devised 
that a proper cut to see a possible signature of the would be to take 1030 MeV for the cut 
and 100 MeV for the Pmin test. The results with these cuts can be seen in Fig. [5l A peak is still 
seen around M^+n ~ 1530 MeV, but of the same size as the statistical fluctuations that one 
observes in other parts of the spectrum. If in an eventual experiment one might decide that a 
cut like the one applied for Fig. |4]is more appropriate to analyze the experiment, one would then 
find the peak around 1530 MeV of Fig. |4]and could see in it the signature of a new resonance. 

After these exercises one still would like to see what happens with large statistics. We have 
made a run with 2000 millions points, out of which around 6420000 points are "good" points 
passing the test of the LEPS cuts. The results are shown in Fig. [B] We see that the narrow 
peak has disappeared in the 1530 MeV region, but a softer bump still remains. A structure of 
this type in an eventual experimental spectrum could easily be interpreted as due to a resonance 
present in the reaction, but in our case it has been created by the cuts taken from the LEPS set 
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Figure 5: M^+„ invariant mass distribution with about 71000 "good" points and with the cuts 
Pmin <100 MeV and M^^ > 1030 MeV. 



up since we do not have any K^n resonance in the model. 

One can visualize the comment made above by separating the mass spectrum of Fig. |6] into 
a smooth background and a "resonant" signal as shown in FiglTl as is usually done in many 
experimental analyses. A structure with resonant shape around 1520 MeV and a width of about 
85 MeV has been generated as a consequence of the cuts. 

In order to see the effects of the cuts we show in Fig. [SI [HI [TOlthe results for the distributions 
which we get by not putting any cut, only the angle cut, or the angle and a 4> cut, i.e., M^+k- > 
1050 MeV, respectively. In figures [51 [5] and [TUj the neutron momentum is the one provided in 
the Monte Carlo runs and is not evaluated according to the MMSA prescription as in the former 
figures. What we observe is that the "resonant like" structure of Figs. [U [7] does not show up in 
any of the runs where the real momenta have been used. Since the (j) cut used in [3] is dependent 
on Pmin and the MMSA prescription is used for the kinematics of the nucleons, we must conclude 
that the combination of these two facts is responsible for the "resonant like" structure seen in 

Figs. E El 



4 Introducing the A(1520) 

Now we introduce in the amplitude — t- K^np a term accounting for the production of the 
A(1520) as explained in section 1^751 The different structure in the momenta of the final states 
for the i4r+A(1520) and 4)p production suggest the incoherent sum of cross sections and we follow 
this option. 
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Figure 6: (Color online) M^+n invariant mass distribution calculated with ~ 6420000 "good" points 
with the same cuts as those made in [3] at LEPS. 



We find that adding the A(1520) production in the amount demanded by the LEPS experi- 
ment does not spoil the mass distribution found in Fig. |4l taking the same seed for the Monte 
Carlo integration. This is shown in Fig. [TTJ However, the introduction of the A(1520) production 
term serves us also another purpose. Indeed, we can make a test to see how the prescription used 
at LEPS to identify resonances works in this case. Since neither the final proton not the neutron 
are detected in the experiment, one must take some prescription to have the best guess for the 
neutron momentum in order to determine the K^n invariant mass of an event. The prescription 
taken at LEPS relies on the concept of minimum momentum as we explained in section \2M\ In 
Fig. [12] we show the results for the K~p invariant mass considering the same cuts as in LEPS. 
We can see that, using the MMSA prescription to calculate the K~p invariant mass as described 
in section \2M one gets a peak at the right mass coming from the A(1520), and the spectrum 
resembles much the spectrum found at LEPS. In this case, if we remove the term of the A(1520) 
production in the — K^K~p amplitude and plot the K~p invariant mass distribution we do 
not find the signal of the resonance, as can be seen in Fig. [T^] (lower area of the spectrum). It is 
interesting to note that, by reasons of symmetry, the spectrum obtained for the K~p invariant 
mass in the absence of any A(1520) term is basically the same as the one obtained for the K^n 
invariant mass in Fig. [SI up to minor differences from the rescattering terms. It is interesting 
to point out here that, should one have a statistical fluctuation around 1530 McV as one has 
in Fig. [31 it would go unnoticed since it would be buried below the A(1520) peak, which has 
become wider than the natural width as a consequence of the MMSA prescription. 

Fig. [12] shows, indeed, that the prescription to identify a resonance using the MMSA pre- 
scription works if the cross section for the production of this resonance is sufficiently large. 
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Figure 7: (Color online) Mj^+^ invariant mass distribution shown in Fig. [U together with the result 
from the sum of two gaussian functions, one peaking at 1520 MeV with a width of 85 MeV and 
another one peaking at 1645 MeV with a width of 185 MeV. 
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Figure 9: (Color online) M^+^ invariant mass distribution calculated with 23900000 "good' 
points with only the angular cut. 
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Figure 10: (Color online) Mj^+„ invariant mass distribution obtained with ~ 12700000 "good' 
events with the (j) cut Mf^+x- > 1050 MeV and the angular cut. 




Figure 11: (Color online) Mj^+^ invariant mass distribution obtained with ~ 66500 "good" points, 
the same cuts as those made in [3J at LEPS and with the A(1520) term. 
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Figure 12: (Color online) Mj^-p invariant mass distributions with and without the A(1520) term 
and making the same cuts that at LEPS. These distributions have been obtained with ~ 6420000 
"good" points. 
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However, let us see what would happen if we actually measured the proton such that we can 
clearly determine the K~p distribution without using the MMSA prescription. The results are 
seen in Fig. [131 We can see that the peak for the A(1520) is sharper than in Fig. [121 We also 
notice that the non A(1520) part of the spectrum (lower area in the figure) changes from one 
figure to the other. This is telling us that the MMSA prescription can indeed be used to identify 
a resonance, but it also distorts the shape of the spectrum and could have some repercussions 
in the K'^n spectrum, as we have already seen. We have checked that the contribution to the 
spectrum at large K~p invariant masses in Fig. [131 comes from production on the neutron, 
where the proton is a spectator. 
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Figure 13: (Color online) M^-p invariant mass distributions with and without the A(1520) term 
and making the angle cut and with M^+j^- > 1050 MeV using the real K~ and p momenta. These 
distributions have been calculated with ~ 12700000 "good" points. 

Before ending the section, we show in Fig. [T31 the equivalent of Fig. [51 using the MMSA 
prescription and the cuts done in [3 , but with the full model, i.e., including (f) and A(1520) 
production. The results resemble those of Fig. [51 but the peak at higher invariant masses is a 
bit more pronounced. This should be compared to Fig. 12 of (3j. 



5 Rescattering contribution 

In Fig. [iniwe show the M^^+n distribution calculated with the MMSA prescription and the same 
cuts as those made in |3l at LEPS obtained considering only the contribution from the single 
scattering, which corresponds to the diagrams Fig. [T^, [TJd, while in Fig. [T51we show the result 
obtained taking into account only the rescattering diagrams of Fig. [1] 
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Figure 14: (Color online) Mj^+n distribution calculated with the same cuts as those done in [3] and 
including the A(1520) term for ~ 6420000 "good" points. 

As one can see, the single scattering terms provide the dominant contribution, although the 
double scattering terms help in increasing a bit the magnitude and their contribution also shows 
some peak around 1520 MeV. 

6 Statistical analysis and comparison with LEPS results 

In this section we address the question of statistics and error analysis. In the first place we 
consider the error from the Monte Carlo simulation. This is given by the standard formula 



where is the relative error for da/dMi„v in a bin of the invariant mass, / is the integrand in 
the integral for the cross section (Eq. (HI])) for the Monte Carlo points that fall inside a certain 
bin, Ne the number of energies of the photon for which we calculate the distributions, rij the 
number of "good" events in one bin for each of these energies and Vj the volume of the phase 
space for each photon energy. 




(70) 
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Figure 15: (Color online) M^+n invariant mass distribution with a statistic similar to the one in 
Fig. SI including only the single scattering contribution and with the same cuts as those done in [3j 



Coming back to Fig. IH we obtain using Eq. (|70|) a relative error of ^ 20 % in the bins of 
invariant mass around 1530 MeV. This is an error that justifies the fluctuation seen in the figure. 

The equivalence of N integration points and N measured events in a experiment from the 
statistical point of view is not straightforward. For this reason, we approach now a different 
procedure which is suited to generate events like in the experiment, which we call simulation 
run. We use the Von Neumann rejection method to generate events. This procedure works like 
a Quantum Mechanical experiment in which events are produced or not according to their prob- 
ability. The Von Neumann rejection method does that: it generates or not an event according 
to its probability measured with respect to the maximum probability, where the probability is 
given by the integrand of Eq. (|44[) . 

We made ten Monte Carlo simulation runs, producing about 2500-3000 events in each of 
them, to simulate the statistics of the 2000 events in ^3j. We can make a statistical analysis of 
these ten different sets of results. The relative error is given by 



at LEPS. 



N' N' 



(^')-^E5f-[(^')-T.9^] 



i=l 1=1 



(71) 



e 



r 




26 







1450 1500 1550 1600 1650 1700 1750 1800 1850 1900 

M^+^ (MeV) 



Figure 16: (Color online) M^+n invariant mass distribution with a statistic similar to the one in 
Fig. m including only the rescattering contribution and with the same cuts as those done in [3j at 
LEPS. 

where N' — 10 and now gi stands for the number of events falling inside a certain bin in each of 
the z = 1, N' samples. 

We show in Figs. [T7ll20l four of the ten sets obtained with the simulation. We also take a 
different bin, 6.28 MeV, to facilitate the comparison with the data of 3 . The error obtained 
with Eq. ([71]) is of the order of 20 % for the bins around Mji+n = 1530 MeV. The runs required 
1210 million Monte Carlo points. We can see that in one of the runs we obtain a statistical 
fluctuation of about 20 counts over the peak at 1530 MeV of Fig. [T31 indicating that, with these 
errors, statistical fluctuations of about the magnitude of the 0+ peak seen in [3] are possible. 
These errors in the simulation are similar to those of Figs. l4l [TT] using the Monte Carlo integral 
method, where similar statistical fluctuations were also observed. It is also interesting to note 
that, like in the integration method, with the simulation we also observe clear peaks around 1530 
MeV in about 1/3 of the runs, and accumulation of strength in that region in practically all of 
them. 

We also can do another test of the statistical significance of our runs. For this we can take 
advantage of the fact that we know the "exact" curve from our large statistic integration, Fig. 
|6] for the </> production model and Fig. [14] for the full model, or the average of the ten Von 
Neumann runs, which is already practically equal to Fig. 1141 Then we can take the standard 
function of any run with respect to the exact curve. The is defined as 
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Figure 17: (Color online) Mj^+^ distribution calculated with the MMSA prescription and the same 
cuts as those done in [3j for 3098 events. 
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Figure 18: (Color online) Mj^+j^ distribution obtained with the MMSA prescription and the same 
cuts as those done in [3j for 2523 events and a different seed than for Fig. [T7]for the simulation run. 
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Figure 19: (Color online) Mj^+„ distribution calculated with the MMSA prescription and the same 
cuts as those done in |3j for 2627 events and a different seed than for Figs. [T71 1181 for the simulation 
run. 
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Figure 20: (Color online) Mj^+„ distribution calculated with the MMSA prescription and the same 
cuts as those done in j3| for 2502 events and a different seed than for Figs. [T71 [TSl [T9] for the 
simulation run. 



31 



Figure 



"exact" reference 



x'/d.o.f 



Fig. a 

Fig. [5] 
Fig. M 

Fig. [m 

Fig. 18] 
Fig. M 

Fig. m 



Fig. E] 

High statistic run with the same cuts 
Fig. M 

Averaged ten Von Neumann runs 



Idem 
Idem 
Idem 



1.50 
2.33 
1.65'= 
0.67 
1.81 
0.56 
1.00 



Table 1: Values of /d.o.f for different figures in the paper. "Three points of the tail of the M^+n 
distribution, around 1800-1850 MeV, of no relevance to the "B"*"" peak have been omitted in this 
evaluation of x^/d-o.f. 



where for a we take ^/Ji- All the curves are normalized to the number of events in the LEPS 

experiment, 1967 events for the K~^n invariant mass. In Table [T] we show the reduced ("we 
have divided the by the number of degrees of freedom (d.o.f) which corresponds to the number 
of bins used) for the different figures that we have used in the discussion. 

We can see that the x^/d.o.f is between 0.56 and 1.81 for the Von Neumann runs. The 
statistical significance of the results of Fig. 21 and particularly. Fig. [TTl which is done with the 
full model, is similar to that of the Von Neumann runs. The /d.o.f for Fig. [5] is a bit bigger 
than in the other cases. This run was obtained using the same seed as in Fig. 2] (equivalent to 
having the same experiment) and making a different cut. This figure is done with production 
alone and not the full model, so the results are only indicative. One can make a run with more 
points to have a smaller x^ /d.o.f but the results do not change qualitatively. However, Fig. [5] 
with the same seed as Fig. |4l a similar number of points, and a different cut is more significative 
because it gives a hint of what might happen in the experiment, having a different shape than 
Fig. |4]and producing larger fluctuations. We do not have the experimental data for this cut and 
hence one cannot compare like in the other cases. 

An interesting test that we can do to compare with the LEPS results is the distribution of 
Pmin- This is done in Fig. [5T] using the full model, i.e, including </> and A(1520) production. 
We can see that we get one sharp peak around zero extending up to about -100 MeV to 100 
MeV. The results are remarkably similar to those obtained in the LEPS experiment (see Fig. 
4 of [3]). Now it is equally interesting to make a plot for the distribution of proton momenta 
in the real case. In Fig. [22] we show the real distribution of the proton momentum component 
along the direction of pmiss. We find a peak similar to the one of Pmin, and a second broad 
peak, around 500 MeV, corresponding to protons which in the CM frame of the np final state 
go in the same direction of Pmiss rather than opposite to it, which is the situation that leads 
to Pmin- In the LEPS analysis these latter events would be given a momentum Pmin and no 
transverse component to the direction of Pmiss- This transverse component is also interesting 
to investigate. In Fig. [53] we show the distribution for the modulus of the transverse proton 
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Figure 21: (Color online) Pmin distribution. 

momentum component, pT- As we can see, the distribution peaks around 50 MeV and stretches 
up to about 200 MeV. 

It is easy to identify the sharp narrow peak of Fig. [22] with events where the proton is 
a spectator (0 production on the neutron), while the broad peak would come from events of 
A(1520) or (j) production on the proton, where the neutron would be a spectator. 

A complement to the further discussion can be seen in Fig. [24l where we show a two 
dimensional plot for the correlation between the K'^n invariant mass evaluated with the MMSA 
prescription and the real one. The events for which the two magnitudes are evaluated have been 
generated using the Von Neumann rejection method. We observe a relatively large dispersion 
of the points with respect to a diagonal line where the two invariant masses would be equal. 
Actually, we see two blocks of points, one around the diagonal and another one to the right 
of the diagonal at large K'^n invariant masses. It is easy to identify the first block with the 
mechanism where the proton is a spectator while the second block corresponds to events where 
the neutron is the spectator. In order to visualize this assertion, we run again the simulation 
removing the A(1520) and (j) production on the proton and the double scattering. We show the 
results in Fig. 1251 As we can see, now there is a very good correlation between the two invariant 
masses and the MMSA prescription should be a good tool to analyze a mechanism of this type. 
Of course, in the actual reaction one has both contributions from 7 scattering with the proton 
and the neutron and the realistic situation of the correlation between the two invariant masses 
is the one given in Fig. [Ml 

The simulation runs, which require a large number of Monte Carlo points to get a few 
thousand of events, are suited to the actual experiment, and one should admit the same errors 
in the experiment of j3] as we have found here. Also one should be ready to admit the possibility 
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Figure 22: (Color online) Distribution for the real proton momentum component along the pr, 
direction. 




100 200 300 

I P^- 1 (MeV) 



400 



500 



Figure 23: (Color online) Distribution for the modulus of the real transverse proton momentum 
component, \pt\- 
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Figure 24: (Color online) calculated using the MMSA prescription versus Mj^+^ obtained 

with the real momentum for the nucleons and the full model, i.e., (j) production on the nucleons and 
A(1520) production on the proton. 
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Figure 25: (Color online) calculated using the MMSA prescription versus Mj^+^ obtained 

with the real momentum for the nucleons omitting the double scattering, (p production on the proton 
and A(1520) production on the proton. 
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to have in the experiment the same statistical fluctuations over the peak of the large statistic 
distribution generated with the MMSA prescription. 

There is one further element in favor of the former argumentation of the statistical fluctu- 
ations. In Fig. 10(a) of [3] we see in the spectrum of the Mj^-p a peak of about 20 events 
per bin around 1650 MeV. This peak is assumed in [3] to be a statistical fluctuation, since in 
Fig. 16 of [3], making a cut M^+k- > 1500 MeV, the peak disappears. Then one should note 
that the peak in M^+„ seen at LEPS around 1526 MeV, when subtracted from the "real" large 
statistics spectrum of our Fig. [Ml corresponds also to about 20 events per bin ( see Fig. 12 
of [3]). The justification of the peak in the M^^p spectrum as a statistical fluctuation has as a 
consequence the acceptance that the peak in the M^+„ can equally be a statistical fluctuation. 
The large statistics calculation done here is important in this respect because it establishes the 
starting point from where to count the fluctuation. Since one already starts from a "resonant 
like" structure around 1530 MeV, the diversion of the data of LEPS over this realistic spectrum 
is smaller than one might assume by inspection of Fig. 12 of ^3 alone. 

7 Conclusions 

We have made a simulation of the "fd — > K~^K~np reaction using a theoretical model that 
accounts for (jj and A(1520) production, the two visible resonance structures in the reaction of 
LEPS [3] . The model takes into account rescattering of one of the kaons with the second nucleon 
in the deuteron. 

A Monte Carlo simulation of the reaction was done and distributions of K'^n and K~p invari- 
ant masses were reconstructed. One of the elements in the LEPS analysis was the consideration 
of a Pmin of the nucleons, by means of which one could minimize the effect of Fermi motion of 
the nucleons in the deuteron and make an educated guess for the final, undetected neutron or 
proton momentum (MMSA prescription). We could see that the procedure allows to recognize a 
resonance if it is present in the theoretical amplitude. However, we also saw that the MMSA pro- 
cedure, together with the cuts to eliminate the cf) contribution, lead to a "resonant like" structure 
around 1530 MeV in the K^n mass distribution, which we could show with a large statistics run. 
Since the peak seen at LEPS is narrower than this "resonant like" structure, we made a thorough 
statistical analysis of the process to see if a statistical fluctuation could be responsible for it. 
For this purpose we made several simulation runs using the Von Neumann rejection method to 
generate events weighted by their probability of production, much as a Quantum Mechanical 
experiment does. The statistical study of these runs with a similar number of events as in the 
experiment [3| provided the intrinsic error. These runs showed that a fluctuation over the peak 
generated by the MMSA prescription and the LEPS cuts around 1530 MeV was possible, with 
a similar size as the peak seen in [5]. On the other hand, in the M^-p spectrum of [3] there is 
a peak around 1650 MeV which was assumed to be a statistical fluctuation since it disappeared 
making a slightly different cf) cut. Then, using the spectrum obtained with high statistics Monte 
Carlo runs we could see that the peak in M/f+„ in over this realistic spectrum has the same 
strength as the peak in M^-p, which was assumed to be a fluctuation in |3]. The parallel study 
done at LEPS of the M^-p spectrum has thus served to set a scale for the statistical fluctua- 
tions in the LEPS experiment, in such a way that one can accept the M^+„ peak as a statistical 
fluctuation like the one admitted in [3] for the peak in the M^-p spectrum which has the same 
strength. 



37 



Another point worth making is that we could see in our Monte Carlo runs that for a same 
seed (equivalent to having a unique experiment) and different cj) cuts, the M^+„ spectra changes 
appreciably, to the point that for some cuts the peak around 1530 MeV looks more prominent 
than for others. This means that, in principle, it is possible to convert an original fluctuation 
into one more pronounced by suitable changes in the cuts, which should be considered as a 
warning when analyzing this type of experiments. In fact, as mentioned above, the use of a ^ 
cut Mj(+K~ > 1050 MeV does not lead to a peak around Mj^-p ~ 1650 MeV in the experiment 
of [3], while the main cut used in [3] gives rise to a peak at this value of the M^-p invariant 
mass. 
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